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A simple model collisionless, dissipative, compressible MHD (Alfvenic) turbulence in a magne- 
tized system is investigated. In contrast to more familiar paradigms of turbulence, dissipation arises 
from Landau damping, enters via nonlinearity, and is distributed over all scales. The theory predicts 
that two different regimes or phases of turbulence are possible, depending on the ratio of steepening 
to damping coefficient (mi/ms). For strong damping (\mi/m2\ < 1), a regime of smooth, hydro- 
dynamic turbulence is predicted. For |rai/m.2| > 1, steady state turbulence does not exist in the 
hydrodynamic limit. Rather, spikey, small scale structure is predicted. 

PACS: 47.65.+a, 47.52.+j 

The theory of compressible MHD (e.g., Alfvenic) turbulence has been a topic of interest for some time §. Alfven 
wave turbulence presents several novel challenges, due to the fact the k-u> selection rules preclude three Alfven-wave 
resonance. Thus, in incompressible MHD, two Alfven waves can interact only with the vortex (i.e., eddy) mode. 
Compressibility relaxes this constraint by allowing interaction with accoustic and ion-ballistic modes (i.e., Landau 
Oh, damping), along with waveform steepening. This naturally leads to the formation of Alfvenic shocklets. Thus, 

S one approach, which is analogous to the noisy-Burgers model in hydrodynamics, is based on the study of nonlinear 
wave evolution equations with external noise drive [e.g., the noisy derivative nonlinear Schrodinger equation (DNLS) 
i equation, in space physics]. Such theories describe turbulence as an ensemble of nonlinear structures, e.g., shocks, 
r— H ' discontinuities, and high-amplitude waves, which are typically observed in compressible (e.g., interplanetary 0]) 
^T 1 , plasmas. This course of investigation was pursued computationally to study the noisy-DNLS equation ||]. Stationarity 
^ was achieved by inserting ad-hoc viscous damping (later linked to finite plasma conductivity jjj) into the otherwise 
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conservative DNLS equation. The DNLS model fails, however, for the important case of j3 ~ 1 (f3 — A-np/B^ is the 



ratio of plasma pressure to magnetic pressure, Bq is an external magnetic field) and the electron-to-ion temperature 



ratio T e /Ti ~ 1 (for instance, in the solar wind plasma), when Alfven waves couple to strongly damped ion acoustic 
modes. As a consequence, the kinctically modified DNLS referred to as the kinetic nonlinear Schrodinger 

equation (KNLS), which exhibits intrinsically dissipative nonlinear coupling, emerges as the superior basic model. 
Numerical solution of the KNLS reveals a new class of dissipative structures, which appear through the balance of 
nonlinear steepening with collisionless nonlinear damping. These structures include arc-polarized and S polarized 
rotational discontinuities J?]], observed in the solar wind plasma and not predicted by other models. The resulting 
• quasi-stationary structures typically have narrow spectra. 

Here, we present the first analytical study of the noisy-KNLS equation as a generic model of collisionless, large- 
CN) ■ amplitude Alfvenic shocklet turbulence. Indeed, this is, to our knowledge, the first structure-based theory of com- 
pressible MHD turbulence in a collisionless system. Stationarity is maintained via the balance of noise and dissipative 
nonlinearity. Dissipation here results from ion Landau damping, which balances the parallel ponderomotive force 
produced by modulations of the compressible Alfven wave train. A one-loop renormalization group (RG) calcula- 



tion (equivalent M to a direct interaction approximation M closure) is utilized. Although the KNLS describes both 



quasi-parallel and oblique waves H, we consider here the simpler case of quasi-parallel propagation. The general case 



will be addressed on future publication. The noisy-KNLS is, thus, a generic model of strong, compressible Alfvenic 
JLT* turbulence and may be relevant to the solar wind, interstellar medium, shock acceleration as well as to compressible 
MHD theory, as a whole. Note that this perspective is analogous to that of the noisy-Burgers equation model of 
compressible fluid turbulence flo|| . Several features which are not common in standard MHD turbulence theories ap- 
pear in this model. It is shown that the dissipative integral coupling renormalizes the wave train velocity, in addition 
to inducing nonlinear damping and dispersion. Moreover, consideration of the resulting solvability condition for a 
stationary state in the hydrodynamic limit (ui, k — ► 0) suggests that KNLS turbulence can exist in one of two different 
states or phases. In the hydrodynamic regime, turbulence consists of large-scale, smooth (lu, k — > 0) waveforms and 
dissipative structures. In the regime when the hydrodynamic limit does not exist, one may expect a small-scale, 
spikey, intermittent (ui,k -/* 0) shocklet turbulence. This hypothesis, however, needs further (e.g., numerical) study. 
The "noisy-KNLS" equation is 



dcfi d<f> d . d\ - 



where <j> = (b x + ib y )/ Bq is the wave magnetic field, / is the random noise, /iq = v A /2Qi is the dispersion coefficient, 
vq is the reference frame velocity, A = 1 is the perturbation parameter, va is the Alfven speed, and fii is the ion 
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gyrofrequency. Unlike the Burgers equation, the KNLS (and DNLS) equation is not Galilean invariant, hence the 
Do term is explicit. The packet velocity vq is renormalized due the broken symmetry between +k and —k harmonics 
induced by Landau damping. This precludes the conventional practice of transforming to the frame co-moving at va 
to eliminate vq. The macroscopic ponderomotive plasma velocity perturbation for a high amplitude Alfven wave is 

- - 1 f°° V 
U 2 = mi \q>\ 2 + m 2 HM% n=- dz', (2) 
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where 7i is the (nonlocal) Hilbcrt operator which represents collisionless (Landau) dissipation. The coefficients mi 
and m 2 are functions of (3 (c s and T e /T,;, only, i.e., 

_1 (l-/3*) + xf(l-/3V7) 
mi 4(l-/3*)2 +x j|(l- / 3V7) 2 ' ( j 

m2 4(l-/3*)2 +x |(l-/3*/7) 2 ' 1 ' 



where (3* = (T e /Ti)/3, 7 = 3 is the polytropic constant, and \\\ = V&P/ttj [(T e /Ti) 3 / 2 exp {(T t - T e )/2Ti}] is the 
parallel heat conduction coefficient which models kinetic collisionless dissipation in fluid models. The term mi 
represents nonlinear steepening of a wave via coupling to the self-generated density perturbation (associated with an 
acoustic mode). The term m 2 corresponds to kinetic damping of a wave by resonant particles, which rapidly sinks 
energy from all harmonics, unlike viscous dissipation. We emphasize that the KNLS is intrinsically a nonlinearly 
dissipative equation, i.e. there is no linear damping retained here. In this regard, we comment that there appear to 
be two meaningful paradigm problems to be explored. One is to retain both linear growth and damping in the KNLS. 
In this case the results will necessarily be quite model dependent. The other, which we pursue here, is to study the 
purely nonlinear problem with noisy drive. This case allows us to isolate and focus on the intrinsically nonlinear 
dynamics of the KNLS equation. 

In Fourier space TC — ik/\k\, so the transformed KNLS is 

(— iuo + ivok + i/iok 2 ) 4> * 

+ i\k*^2 [ <f> k' <t> k" i> k-k'-k" [mi + im 2 sign(fc - k')} J = fk , (4) 

where the function sign(a;) = x/|a;|. The stochastic noise fk is assumed to be zero-mean, and 5-correlated in space 
and time. To extract information from Eq. (0), we utilize the direct interaction approximation (DIA) closure |p|, pd] , ^2[ . 
We follow the approach used in jl2| which is different from (llj , where the correlation functions of cubic Schrodinger 
turbulence were found. We emphasize the statistical nature of our analysis. Indeed, turbulence renormalizes the 
coefficients of the evolution equation, leaving its functional form unchanged. The renormalized phase velocity v and 
dispersion /1 (see discussion below) appear in the second order in A. Since the KNLS (and DNLS) is not Galilean 
invariant, the vertex A is also renormalized. This is a third order effect. Because of mathematical difficulties, we 
exactly calculate v and fi, only, and for A, we provide a simple heuristic argument. In general, the noise finction / is 
also renormalized at higher orders in A. Usually, in one-loop RG analyses, such turbulent corrections are assumed to 
be small. Note that knowledge of turbulent A and / is not required for predicting the existence of the hydrodynamic 
limit (see below). 

We expand <j)k in a power series with respect to the perturbation parameter A: <pk — <pY + A^^ + A 2 0^ + ••• 
and equate terms, order by order, in A. To second order, we have 
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where the bare propagator Gq(uj, fc) = i/(u> — kvo — fc 2 Mo)- ln the DIA, we take <ffy ~ <f>$ '. The terms proportional 

to k and k 2 in the left hand side act to modify vq and /iq. Thus, the nonlinear term of Eq. (||) represents an 
amplitude dependent correction to both the velocity and dispersion coefficients, and Eq. (|E|) is a recursive equation 
for the renormalized coefficients v and fi. The fixed point of this recursion relation gives the self-consistent values of 
these coefficients. Replacing the bare vq, /xo with their amplitude dependent counterparts v, fj., we write 
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We should note that v and /x will now assume complex values, 



IVi, 



fi = fx r + ifii. 



(7) 



The real parts, v r and fj, r , represent the amplitude dependent speed of a wave packet (note, there is a momentum 
transfer from waves to resonant particles in this model) and nonlinear dispersion (i.e. an amplitude dependent 
frequency shift in Fourier space), respectively. The imaginary parts, Vi and /ij, correspond to damping processes. In 
particular, Vi describes the exponential damping (a la phase mixing of a wave packet) and fii describes turbulent, 
viscous dissipation. It is easily seen that for W2 — > the KNLS may be written in the co-moving frame with 
vq = 0. Thus, no additional phase-mixing terms appear, since the terms k'k and k"k vanish upon integration over 
-co < {k',k"} < oo, in the hydrodynamic limit. The collisionless damping breaks this symmetry of the +k and — k 
parts of the spectrum, thus resulting in the novel phase-mixing and phase velocity renormalization terms (analogous 
to nonlinear frequency shifts) encountered here. 

We seek solutions in the hydrodynamic limit to — > 0, k — > 0. For simplicity, we assume for noise the white 
noise statistics, i.e., /* = /. This assumption is not too artificial, since MHD waves are usually pumped at large 
scales (small-fc) and the large-fc tail is heavily damped by collisionless dissipation, which is an increasing function 
of k. Ordered by powers of k, the nonlinear term in the integrals contains fc 3 , fc 4 , . . . contributions. However, the 
hydrodynamic behavior is completely determined by the small-fc, u> limit. Thus, by omitting higher-fc terms, Eq. (|^) 
naturally splits into two equation for v and respectively. The uj' , ^''-integrations can be easily performed in complex 
plane. It is convenient to introduce dimensionless variables x' — k'fi/v, x" — k" ' ji/v. The fc', fc"-integrals in Eq. (|^), 
diverge as fc', fc" — > (i.e. infrared divergence). The integrations can be performed consistently only in the limit 
where the infrared cut-offs satisfy the inequality x' c ,x" <C 1. Quite lengthy, but straightforward complex integrations 
yield 



v r + IVi 



A 2m m 



(2tt) 2 v 3 



[i r + ifn =-f 



(2tt) 2 v 4 



In x c + mln(l - 2p) + -j- I — + '! 



(8a) 



4 A 2 m/j, | „ lna; r . 



mw(/i — 3) + mv{h — 3/i) 



In 2 x c m 1 + 2v(l — /i) +lna; c 2m{\x — v) 



m(4/i-5u + l) +F(v,(m) 



(8b) 



where v = v — v* , jl = fi — /i*, m = mi + im2, m — [m — m*] sign(u^i), m = m + m* , v — v/v, /2 = fJ,/ fl, and the 
dimensionless infrared cut-off is x c = k m i n fii/vi. The function F(v,p,) is positive definite and contains no explicit 
divergences ~ x c . Since we are concerned with the hydrodynamic limit, where i c < 1, the detailed structure of this 
function is not significant. Returning to standard notation, extracting real and imaginary parts, and keeping the 
leading, divergent (in x c ) terms, we have from Eq. ( j3a| ) 
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where a = sign(/z r //ii). As can be easily seen, v r /vi ~ l/(ln 2 k m i n ) — > as k m i n — > 0, so we obtain from Eq. (8b) 

l2 ,„ / //. 
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For the coefficients w and we may now write 
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Note that the factor lna; c = \n(k m i n fj,i/vi) ~ ln(lrifc m j n ) makes an insignificant cut-off correction. Nonzero v r arises 
due to wave momentum loss via interaction with resonant particles and reflects the process whereby a nonlinear 
wave accelerates in the direction of steepening (i. e. v r > for \3 < l,T e = Tj), an effect which is observed in 
numerical solutions of the KNLS equation JtJ. This effect is logarithmic for k m i n — > 0. Negative Vi corresponds 
to exponential damping due to phase mixing, and is proportional to the dissipation rate m 2 . The coefficient /x r 
represents turbulent dispersion, and the coefficient fix < corresponds to turbulent viscous damping. By analogy 
with noisy-Burgers equation |]l2| , Eqs. ( |l0| ) for the turbulent transport coefficients yield the pulse propagation scaling 
exponents for the hydrodynamic regime, which are defined by divergences at the cut-off. For diffusion term, we have 
5x 2 /St ~ ^ ~ \Sx\ (~ 1/fcmm), so that \Sx\ r~j St. This corresponds to symmetric ballistic dispersion of the shocklet 
waveform. For the velocity term, we write (as v r — ► when fc m i„ — ► 0) x/t ~ Wj ~ VAlnx c ~ const, that is a; ~ i. 
This corresponds to ballistic translation of the shocklet. 

We now construct the quantity x c = k m i n fii/vi from Eqs. (||) to determine when our cut-off approximation i c <Cl 
is valid. Note that i c <1 must be satisfied for a self-consistent, hydrodynamic regime solution. Dividing Eqs. (9b,[k]) 
by Eq. (^), we derive a system of equations which is easily simplified to give the condition 
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< 1, i.e. 



Again, for x c <C 1, we may omit the small x c term. This equation has maximum at (x c ) max = 
solution of this equation for small x c exists only when x c < (x c ) max . When x c > (x c ) max , Eq. (11) does not have 
a small- a; c solution, so no stationary state is possible in the hydrodynamic limit. To clarify the physical meaning of 
the control parameter x c , we write it as x c = (k^nHi) I {kminVi) ■ Obviously, x c is just a measure of the efficiency 
of turbulent viscous damping ~ k 2 relative to collisionless (Landau) damping (distributed in all scales). Smallncss 
of x c indicates a situation of stronger Landau damping and weaker linear turbulent (viscous) dissipation. The two 
cases of x c lesser or greater (x c )max thus correspond to different states of turbulence. The regime of hydrodynamic 
turbulence [i. e. x c < (x c ) max ] corresponds to strong damping, |m 2 | 3> \ttii\, which dominates nonlinear steepening. 
Large-scale waveform structures are possible, consistent with the notion of a hydrodynamic regime. The turbulent 
viscous damping dominated nonlinear dispersion in this case, fir/ n i — 4mi/3m 2 -C 1. The opposite regime of 
"shock" turbulence (i. e. x c > (x c ) max ) corresponds to weakly damped Alfven waves, |m 2 | -C |mi| (however, Landau 
damping still dominates the small-scale dissipation), where a stationary, hydrodynamic regime is not possible. In this 
case, nonlinear steepening is balanced by turbulent dispersion, resulting in a state of small-scale coherent nonlinear 



structures, steep fronts and discontinuities. The bifurcation point can easily be found from [x c ) r 
(0) as 
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(the exact numerical solution yields ~ 1.3). The coefficients mi,m 2 depend on plasma parameters, i. e. on p and 
T e /T t . We plot the condition Eq. @ in the form of a vs. T e /Xi diagram m Fig. |. The region inside the curve 
corresponds to |mi/m 2 | < 1, i.e. a phase of hydrodynamic turbulence. The outer region corresponds to a phase of 
bursty turbulence of steep nonlinear Alfven waves. 

For completeness, the perturbation parameter A must be renormalizcd, because the KNLS is not Galilean invariant. 
Corrections to A follow from the third order expansion and so laborious that they are left for a future publication. 
We can, however, estimate the renormalized A as follows. The energy spectrum is 
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The fluctuation level 
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should be independent of the cut-off x c , thus A ~ In" 1 cc c . Of course, the fluctuation level may only depend on the 
noise strength, /, and the dissipation rate, m 2 . As is expected, B 2 varies as m^ 1 while noise is constant. 

To conclude, we have presented the first analytical theory analysis of a noisy-KNLS (and DNLS) model. The 
noisy-KNLS describes turbulence of kinetically damped (at (3 ~ 1) nonlinear Alfvcn wave turbulence, i.e. a turbu- 
lence of dissipative structures 0, discontinuities and shock waves. The renormalized wave velocity and dispersion 
coefficients, as well as the pulse propagation exponents, were calculated. Two different phases of turbulence were 
identified, depending on the nonlinearity-to-dissipation coefficient ratio, m\jmi. For \m\lm2\ < 1 a stationary state 
of hydrodynamic (k — > 0, w — > 0) turbulence (with noise) is predicted, while for \mi/m^\ > 1 such a state is precluded 
and small-scale bursty, spikey turbulence is indicated. A phase diagram in the space of f3 and T e /Ti is given. These 
findings may be pertinent to recent observations of multiple states in solar wind plasma turbulence. 
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FIG. 1. The /3-Te/Ti-diagram of state. The region inside the curve corresponds to highly damped turbulence. No steep 
fronts appear. There is wave steepening in the region outside the curve. 
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